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ARCSECOND SCALE KINEMATIC AND CHEMICAL COMPLEXITY IN CEPHEUS A-EAST 
C. L. Brogan', C. J. Chandler', T. R. Hunter', Y. L. Shirley-, A. P. Sarma^ 
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ABSTRACT 

We present results from SMA observations of the star forming region Cepheus A-East at ^ 340 GHz 
(875 fim) with 0."7-2" resolution. At least four compact submm continuum sources have been detected, 
as well as a rich forest of hot core line emission. Two kinematically, chemically, and thermally distinct regions 
of molecular emission are present in the vicinity of the HW2 thermal jet, both spatially distinct from the submm 
counterpart to HW2. We propose that this emission is indicative of multiple protostars rather than a massive 
disk as reported by Patel et al. (2005). 

Subject headings: ISM: individual (Cep A) — ISM: lines and bands — ISM: molecules 



1. INTRODUCTION 

The Cepheus A-East (hereafter CepA) star forming region 
lies at a distance D 725 pc and has Lboi ^ 2.2 x 10"* L0, con- 
sistent with a cluster of BO. 5 or later stars (Blaauw, Hiltner, 
& Johnson 1959; Sargent 1979; Mueller et al. 2002). At cm- 
wavelengths, CepA consists of several compact sources called 
HWl, HW2, ... HW9, which He along a roughly inverted Y- 
like structure (e.g., Hughes & Wouterloot 1984; Garay et al. 
1996). It is currently unclear how many of these compact ion- 
ized structures correspond to individual protostars. For ex- 
ample, much of the cm-wavelength emission from the HW2 
region is due to a bipolar thermal jet rather than a Stromgren 
sphere (Rodriguez et al. 1994; Curiel et al. 2006). 

A wide range of other signposts of on-going star forma- 
tion have been observed in CepA including several outflow 
components (e.g. Codella et al. 2005, and references therein), 
however the locations of the powering sources remain uncer- 
tain. On smaller sizescales copious OH, H2O, and CH3OH 
maser emission has been detected toward several of the cm- 
wavelength sources (Vlemmings et al. 2006, and references 
therein). Patel et al. (2005) report the detection of a mas- 
sive molecular gas and dust disk toward the HW2 source from 
Submillimeter Array'* (SMA) observations of methyl cyanide. 

In this Letter, we present new and archival SMA 345 GHz 
observations toward CepA (including those of Patel et al. 
2005), concentrating on the continuum data and the spa- 
tial, kinematic, and temperature information provided by the 
wealth of spectral lines; complete details of the spectral line 
results will be presented in a future paper. Our observations 
and results are presented in § 2 and are discussed in § 3. 



TABLE 1 
SMA Observing Parameters 



Date 


u-v range 


tint 


USB/LSB° 


Line Beam* 




(kA) 


(hours) 


(GHz) 


(" X " (°)) 


30 Aug 2004'' 


20-190 


5.7 


331.4/- 


0.9 X 0.8 (77) 


26 Sep, 18 Oct 2004 


14-130 


5.0 


343.0/333.0 


1.9 X 1.2 (71) 


05, 07 Oct 2005 


10-80 


4.1 


346.6/336.6 


2.0 X 1.9 (33) 



"Approximate center frequency of 2 GHz wide sidebands. 

''Synthesized beam and position angle of the USB line data. 

■^Archival Patel et al. (2005) data, only the upper sideband was analyzed. 



package. The quasars BL Lac and 3C 454.3 were used for 
phase calibration. 3C 454.3, Uranus, and Saturn were used 
for bandpass calibration. Comparison of the amplitude cal- 
ibration (from the quasars) applied to Uranus vs. a model 
of its baseline dependent flux density suggests that the CepA 
amplitude calibration is accurate to within ^15%. 

Continuum subtraction, imaging, deconvolution, and self- 
calibration were carried out in AIPS. After extracting the con- 
tinuum using line-free channels in the u-v data, the continuum 
for each observation/sideband were separately self-calibrated; 
the derived phase and amplitude corrections were also applied 
to the line data sets. The spectral line data were Hanning 
smoothed during imaging to produce a final spectral resolu- 
tion of ^ 1.4 km s"' (except for the archival SMA dataset 
with 5.9 km s"^ spectral resolution). To create the highest 
possible sensitivity continuum image, all of the final contin- 
uum datasets were combined in the u-v plane and imaged. 



2. OBSERVATIONS & RESULTS 

The SMA observing parameters are provided in Table 1. 
The data were taken with a channel width of 0.8125 MHz 
(^ 0.7 km s"') except for the archival data which have 8 
times poorer spectral resolution across most of the band 
(the CH3CN (18-17) K=0-3 fines were observed with - 0.7 
km s"' resolution). All five epochs employed seven anten- 
nas. The data were calibrated using the MIRJAD software 

' National Radio Astronomy Observatory; cbrogan@nrao.edu 
^ University of Arizona; Bart J. Bok Fellow 
' Physics Department, DePaul University 

■* The Submillimeter Amy is a joint project between the Smithsonian As- 
trophysical Observatory and the Academia Sinica Institute of Astronomy and 
Astrophysics. 



2.1. CepA Submillimeter Continuum 

Fig. la shows the combined naturally weighted 875 /im 
SMA continuum image with a beam of l."3 x l."l and inte- 
grated flux density of 7.7 Jy; the rms noise is 13 mJy beam"' . 
The peak of the single dish submm source (see e.g. Mueller 
et al. 2002) has been resolved for the first time into a num- 
ber of distinct components. Archival 2002 Very Large Array^ 
(VLA) 3.6cm continuum contours are also shown on for com- 
parison. While there is rough agreement between the peak 
875 fjm continuum emission and the well-known ionized ther- 
mal jet HW2 (e.g. Rodriguez et al. 1994; Garay et al. 1996), 

^ The National Radio Astronomy Observatory operates the Very Large 
Array and is a facility of the National Science Foundation operated under a 
cooperative agreement by Associated Universities, Inc. 
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Fig. 1 . — (a) Combined SMA 875 fim continuum image (greyscale and black contours) witli a resolution of lf'3 X 1"0 (P.A.=79°) and contour levels of -40, 
40 (3a-), 80, 120, 200, 300, 400, 600, 1000, and 1600 mJy beam"' . This image has not been coirected for primary beam attenuation. Green VLA 3.6 cm contours 
at 0.06 (3cr), 0.12, 0.2, 0.3, 0.4, 0.5, and 1.5 mJy beam"' and 0"23 resolution are superposed. The region shown in Fig. lb is indicated by the dashed box. (b) 
Superuniform weighted 875 fim contour map (black) created using baselines longer than 40 kA and restored with a Of'6 beam; the contour levels are -50, 50 
(4a), 100, 150, 300, 500, 700, and 900 mJy beam"'. The green contours are the same as in Fig. la. The synthesized beams are shown in the lower left. Colored 
crosses mark the peak positions of the molecular species listed. Prominent cm-A and submm sources are also labeled. Sample spectral line profiles from the 
positions indicated on Fig. lb for (c) SMA2 and (d) HW2-NE. The displayed transitions are C^*S (7-6), CH3OH (14, 8/l±-156,9A±), NHoCHO (I62J5-I52.14), 
H2CS (10o.io-9o.c)), HC3N (38-37), SO2 (I67.9-I76 12), and SO (llio-lOio). Note that at the spectral line angular resolution of these transitions (~ 2") the two 
positions are not completely independent. 



there is excess emission to the NW of HW2 suggestive of ad- 
ditional unresolved structure. 

In order to better compare the morphology of the cm 
and submm emission, we have also created a superuniform 
weighted image using baselines > 40 kA, resulting in a beam 
size of 0"8 X 0"7. To further delineate the morphology of 
the emission in the vicinity of HW2, in Fig. lb we present 
the superuniform weighted image restored with a 0."6 beam 
(equivalent to 1.8 times the longest baseline sampled) which 
emphasizes the locations of the clean components. Fig. lb 
reveals the presence of at least two distinct submm sources 
in the vicinity of HW2, HW2-SMA and SMAl, in addition 
to an extension NW of HW2-SMA which we denote SMA2 
(see §3). A weak extension south of HW2-SMA (denoted 
SMA3) is also visible, which may be a submm counterpart 
to the low mass protostar VLA-R5 (Curiel et al. 2002). The 
combined morphology of HW2-SMA, SMA2, and SMA3 is 
the structure reported by Patel et al. (2005) as a dust disk; their 
Fig. 1 with ^ 0."75 resolution also shows that the centroid of 
submm emission is not centered on the axis of the HW2 jet. 
The HW2-SMA 875 ^m continuum peak is within 0."1 of the 
proposed location of the powering source of the HW2 thermal 
jet based on high resolution (0."05) 7 mm VLA data (Curiel 
et al. 2006) which is well within our absolution position un- 
certainty of 0."15. Two additional submm cores are detected 
to the south of HW2, one coincident with the cm- wavelength 
source HW3c (denoted HW3c-SMA) and another located at 
the NE tip of the cm-wavelength source HW3b (designated 
SMA4). No distinct compact 875 fim counterparts to HW8, 
HW9, HW3a, HW3b or HW3d are detected in the SMA data. 

2.2. CepA Hot Core Line Emission 

Within the 8 GHz of total bandwidth observed (Table 1) 
we detect more than 20 distinct species along with a num- 
ber of their isotopologues. A few transitions with low ex- 
citation energies CS(7-6), \\^^C0^(A-3)) show 
extended emission over the full primary beam (^ 36") corre- 
sponding to the previously studied outflows (see e.g. Codella 
et al. 2005). However, the emission from most species is 
restricted to compact regions (< 2") that coincide with the 



submm continuum emission (Fig. lb). At l"-2" resolu- 
tion (Table 1) the species exhibiting compact emission in 
the vicinity of HW2 are strongest at one of two distinct ve- 
locities: -5.0 ±0.5 or -10.5 ±0.5 km s"' (Fig. lb,c,d; see 
also Codella et al. 2006). These two kinematic features are 
also spatially distinct. Without exception, molecules that are 
strongest at ~ -10.5 km s"' have peak positions ^ 0"25 (200 
AU) E/NE of HW2-SMA (we denote this position HW2-NE), 
and molecules that are strongest at ^ -5 km s"' peak at the 
position of SMA2. Although this dichotomy exists for all 
observed species, a few abundant high density tracers like 
CH3CN and C^'^S, show emission of nearly equal strength to- 
ward both positions (Fig. lc,d). 

Using 3mm Plateau de Bure (PdBI) data Martm-Pintado et 
al. (2005) also find that SO2 emission peaks to the E/NE of 
HW2 at a velocity of ^ -10.5 km s~' and suggest that this 
emission is due to a distinct intermediate mass protostar; the 
SMA and PdBI SO2 positions agree to within the absolute po- 
sition uncertainty of 0"15. This result has recently been con- 
firmed by the VLA detection of SO2 emission and weak 7mm 
continuum emission at the position of HW2-NE (Jimenez- 
Serra et al. 2007). In an extensive PdBI spectral line study 
Comito et al. (in preparation) also find strong spatial, chem- 
ical, and kinematic differentiation in general agreement with 
the SMA results. The two velocity components at ^ -5 km s~' 
and ^ -10.5 km s"' have also been observed in single dish 
H2CS, CH3OH, and HDO data with resolutions ranging from 
10" to 30" (Codefla et al. 2006). The single dish emission 
from both velocity components is extended. Our SMA data 
are insensitive to spatial structures > 15", but it seems likely 
that SMA2 is associated with the larger scale ^ -5 km s"' 
component. 

For CH3OH towards SMA2 (at - -5 km s"'), and SO2 
and HC3N towai-d HW2-NE (at - -10.5 km s"'), we have 
measured enough transitions to construct rotation diagrams 
(Fig. 2; see Goldsmith & Langer 1999). Emission from 
the two velocity components are well-separated, even though 
they are not completely spatially resolved by our observations 
(Fig. lc,d). The CH3OH and SO2 data have been corrected 
for optical depth effects by first estimating the opacity of 
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TABLE 2 

Derived Parameters of Submillimeter Cores 



200 400 600 800 1000 

E„ / k (K) 

Fig. 2. — Rotation diagrams for CH3OH A (•) and E (o) transitions toward 
SMA2 and for SO2 (□) and HC3N (A) toward HW2-NE (see Fig. lb). The 
CH3OH and SO2 data have been corrected for optical depth effects. 

one transition using a less abundant isotopologue and assum- 
ing Galactic abundance ratios of ^^C/^'^C=70 and 32s/34s=23 
(Milam et al. 2005; Chin et al. 1996, and references therein). 
Then, assuming LTE, the other transitions were corrected us- 
ing the formalism described in Sutton et al. (2004). The 
CH3OH optical depths are significant (up to 28.2); the SO2 
optical depth is more moderate (up to 2.4). A moderate optical 
depth in the lower lying HC3N Unes is also possible (see e.g. 
Wyrowski et al. 1999), correction for which would yield an 
HC3N Trot more consistent with SO2. Three of the lower en- 
ergy A-type transitions of CH3OH and the two lowest energy 
transitions of HC3N are also detected toward HW3c-SMA and 
SMA4; rotation diagram analysis for these two submm cores 
yield 7;„,=65±25 K. 

In addition to the kinematic and chemical dichotomy be- 
tween SMA2 and HW2-NE, their T,-„, differ by more than 100 
K. Both values of Trot (Fig- 2) are significantly larger than 
reported by Torrelles et al. (1999) based on ~ 1" resolution 
VLA NH3 data (r,.„,=30-50 K), influenced by their reported 
non-detection of NH3 (4,4) [E,/k = 202 K]. An NH3 (4,4) in- 
tegrated intensity image from our reduction of these archival 
data is shown in Figure 3; emission is clearly detected toward 
both HW2-NE and SMA2. A new analysis of the NH3 data 
gives Trot within a factor of two of those reported here. Our 
Trot are also larger than those estimated by Patel et al. (2005) 
(25-75 K) though the sense of the temperature gradient is in 
agreement (i.e. warmer in the E than W). Although Patel et al. 
(2005) report that CH3CN (18-17) is not detected above K=3, 
our reduction of their archival SMA data detects CH3CN up 
through K=8 (Eu/k = 607 K, K=6 is shown in Fig. 3), and we 
find similar r„M to those shown in Fig. 2. Using 30m data, 
Martin-Pintado et al. (2005) derived T^c- 150 K for 3mm 
SO2 and HC3N transitions toward HW2-NE, however, beam 
dilution may well play a role in this single dish result. 

3. DISCUSSION 

The observed distribution of submm continuum emission in 
the innmediate vicinity of HW2 (±0.5") allows for two pos- 
sible interpretations: a single elongated structure or two (or 
more) marginally-resolved individual sources. Although the 
kinematic dichotomy between the positions labeled HW2-NE 



Name 


R.A., Dec." 


St 


Tb 




Mgas 




(s,") 


(Jybeam ') 


(K) 


(K) 




SMAl 


17.873, 50.56 


0.28 


5.0 


20-100 


1.0-0.1 


SMA2 


17.927, 49.74 


0.80 


14.4 


115-135 


0.32-0.27 


HW2-SMA 


17.999,49.43 


1.17 


21.0 


40-100 


2.0-0.6 


HW3C-SMA 


17.949,46.07 


0.25 


4.5 


40-90 


0.3-0.1 


SMA4 


17.820, 45.33 


0.16 


2.9 


40-90 


0.2-0.1 



^J2000 coordinates added to 22h56m, 62°01', the relative and absolute po- 
sition uncertainties are better than ()"02 and 0" 15, respectively. 
''Using superuniform weighted image with restoring beam 0"84 X Of 70. 



and SMA2 (Fig. lb) could be interpreted as a velocity gradi- 
ent across a continuous structure (e.g. Patel et al. 2005), the 
dramatic chemical and thermal differentiation demonstrated 
by our multi-species analysis is difficult to explain with this 
picture. Chemical and thermal gradients might be expected 
in the radial and vertical directions within a protostellar disk, 
but such dramatic azimuthal asynmietries are more difficult to 
conceive in a single source scenario. We also note that beyond 
a velocity gradient, the kinematic evidence for a Keplerian ro- 
tating disk is quite weak. For example, the weakness of the 
emission at the central HW2 position in position- velocity (P- 
V) diagrams such as those in Patel et al. (2005, and in our own 
data, but not shown) are inconsistent with theoretical expecta- 
tions unless a central hole is invoked (e.g. Richer & Padman 
1991). Even with this modification, it remains difficult to ex- 
plain the unequal position offsets of the two velocity peaks 
from the HW2 stellar position (see NH3 in Fig. 3 and the K=3 
P-V diagram of Patel et al. 2005). In contrast, the presence 
of multiple sources at different velocities naturally explains 
the observed behavior. We therefore favor the multiple source 
hypothesis (see also Martin-Pintado et al. 2005), with at least 
fliree sources in flie vicinity of HW2 (HW2-SMA, HW2-NE, 
and SMA2). The remaining discussion proceeds with this in- 
terpretation. 

Table 2 summarizes the properties of the submm cores iden- 
tified in Fig. lb (excluding SMA3). The gas masses were es- 
timated using 



RS,yD Tdusi 



gas ■ 



B[v,Td\ (l-exp[-Tto]) 



(1) 



and assuming a gas-to-dust ratio R= 100, D = 725 pc, and 
a dust opacity nm^m = 1-84 cm^ g"' extrapolated from Os- 
senkopf & Henning (1994) for thin ice mantles and density 
10^ cm~^. Values for the peak flux density Si, (none of the 
cores appear to be resolved), the continuum brightness tem- 
perature (Th), and the range of assumed dust temperatures (7^) 
are also listed in Table 2. The Mgas have been corrected for the 
continuum opacity calculated from rj„5, = -ln[l -(7/,/r^)]. 

The derived masses are very sensitive to the assumed tem- 
peratures (Table 2). For SMA2, HW3c, and SMA4, we have 
used the range of dust temperatures from rotation diagram 
analysis (§2.2). Temperafiires for HW2-SMA and SMAl are 
difficult to estimate due to the proximity of strong emission 
from HW2-NE and SMA2. No species peak at the posi- 
tion of HW2-SMA suggesting it is not very warm (i.e. com- 
pared to SMA2 or HW2-NE), though its association with the 
strong bipolar jet suggests it is unlikely to be very cold ei- 
ther. Thus we have assumed a moderate temperature range of 
40 to 100 K for HW2-SMA. Since SMAl is lacking both a 
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Fig. 3. — Integrated intensity images (over tlie full width of the observed emission) of archival VLA NH3 (4,4) data with If 1 X If'O resolution and archival 
SMA data with Of'87 X 0"75 resolution. White 875 ^m contours from Fig. lb, and a green 0.6 mJy beam"' 3.6cm contour are superposed. The black + symbol 
shows the proposed location of the star powering the HW2 thermal jet (Curiel et al. 2006). The beams are shown in the lower left and the reference position is 
22^56™ 17.998", 62°01'49"50 (J2000). 



cm-wavelength counterpart and any distinct line emission we 
assume a cooler lower limit (20 K). With these assumptions 
the total gas masses are fairly low, ranging from 0.1 to 2.0 
Mq. If these are preprotostellar cores, this result implies that 
a massive star will not form; if, however, an embedded star 
is present (almost certainly true for HW2-SMA) these masses 
are consistent with typical values observed for intermediate 
mass protostellar disks (e.g. Hamidouche et al. 2006). 

No compact submm emission corresponding to HW8, 
HW9, HW3a, HW3b, or HW3d is detected. HW3a, HW8 and 
HW9 are variable at cm-wavelengths and are thought to be 
low mass pre-main-sequence stars (e.g., Hughes 1988; Garay 
et al. 1996), hence it is unsurprising that any submm emission 
is below our 3 a detection threshold of 40 mJy beam"'. Based 
in part on the detection of very compact (0."1) 2 cm emission 
towai-d HW3d (Hughes 1988), Garay et al. (1996) suggest that 
HW3d contains its own internal energy source but also has 
thermal jet-like extended cm-wavelength emission. Given its 
lack of distinct submm continuum or line emission we sug- 
gest instead that HW3d is a one-sided ionized jet emanating 



from HW3C-SMA (see Fig. lb). HW3b is consistent with 
being a one-sided jet emanating from SMA4 or possibly the 
counterjet to HW3d (also see Garay et al. 1996). 

We detect at least five submm sources (HW2-SMA, SMAl, 
SMA2, HW3C-SMA, and SMA4) within a projected radius 
of 4" (2900 AU). If the five low mass sources HW3a, HW8, 
HW9, VLA-R5, and VLA-R4 (Hughes 1988; Curiel et al. 
2002) are included and equal clustering in the perpendicu- 
lar dimension is assumed, then the implied protostellar den- 
sity is 5.7 X 10^ pc""*. This approaches the minimum theoreti- 
cal value (10^ pc~^) needed to test the induced binary merger 
hypothesis proposed as a formation mechanism for the most 
massive stars (Bonnell & Bate 2005). 
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by start-up funds to A.P.S. at DePaul Univ. 
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